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The development of structured methods for proportional-integral (PI) 
controller design for systems with time delay are proposed in this article. 
Several PI controller design methods for time-delay systems have been 
reported. However, combining two or more methods to form new ones have 
not been given serious attention. The system stability region in the controller 
parameters space was determined by plotting the stability boundaries. In this 
study, the controller gains were first obtained using genetic algorithm (GA), 
weighted geometric center (WGC), and centroid of convex stability region 
(CCSR). Thereafter, these gains were combined by finding the centroids of 
lines joining any of the two gain locations, and triangle whose vertices are 
the location of the three gains in the convex stability region, thus yielding 
four additional methods, M1, M2, M3, and M4. Compared to a particular 
existing method, some of the proposed methods yield faster response speed 
at the expense of reference input tracking, while the reverse is the case for 
others. Any of the proposed methods (M1, M2, M3, and M4) can be selected 
depending on the system performance specifications. 
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1. INTRODUCTION 

Many practical control systems are modelled as time-delay system because of the negative effects of 
time delay on their normal operations [1]-[3]. Time delay can arise from signal processing overhead incurred 
by the system components and signal transmissions in the control loop [4], [5]. Therefore, designing a 
controller for this class of system requires full knowledge of how the system performance can be adversely 
influenced by time delay, and even be distabilised in worst case [6]. Among the innovations that contributes 
to the field of instrumentation and control in the 20" century, proportional integral derivative (PID) controller 
was rated second [7]. PID controllers are commonly used in industrial applications [8] due to their popularity 
[9], simple algorithm [10]-[13], robustness [12], [14], [15], simple design methodology and effectiveness in 
system control [16]. Most a time PID controller is implemented in the form of proportional-integral (PD 
because its derivative component is not normally use [17]. 

The first step in the PI controller design for systems with time delay is establishing the system’s 
stability boundary in the space of controller parameters [18], being the most common modern method in this 
type of controller design [19]. To design PI controller for time-delay system, [19] proposed a method that 
was based on plotting the stability boundary locus in controller parameters plane and then use the centroid of 
convex stability region (CCSR) method for determining the appropriate controller gains. The superiority of 
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this method over the existing methods was established via simulation and experiments. Consequently, 
analytical method for determining the CCSR in the parameter space of PI controller in the control of 
integrating systems with time delay [17] and unstable or stable first-order plus dead-time systems [20] were 
proposed. However, because of analytical approximation involved, the method sacrifices accuracy for non- 
repeated plotting of stability boundary locus each time new system is considered. In kp-k; stability region for 
a fixed kg, (where kp, k; and kg are proportional, integral and derivative gain respectively) was actually 
established for designing of PID controller, but later adopted for PI controller by allowing kg = 0, to yield 
stability region in k,-k, stability region [21], [22]. Using this region and genetic algorithm (GA), optimum 
kp and k; were determined. A PI controller design method, which involves determination of system stability 
region in the controller parameters space and the centroid of the region was developed by [23]. The required 
controller is the coordinate of the centroid of the determined region. The method is simple, gives good time 
domain performance, and saves implementation time, compared with some common methods. 

In the design of PI-PD for time-delay systems by [24], the inner PD loop was designed by 
determining the k,-kg stability region, weighted geometric center (WGC) method was used to find the controller 
gains within the region and based on these gains, the loop was reduced to a single block. After the reduction 
of the inner loop, the PI was then designed in the k,-k; stability region via WGC. The method is simple, 
robust and gives good performance compared to available methods in the literature. Onat [25] proposed a 
method of PI design that involved determination of k,,-k; stability region and then calculated the WGC of the 
region whose coordinate values were the required gains. From the simulation results, the method yielded 
good compromise among the TDPMs. A method of PI controller design which rely on generating the system 
stability region in the k,-k; plane and then calculate the region WGC was present by [23]. This method is 
simple and reliable which make it to be very useful for control system designers. 

Several methods of PI controller design for time-delay systems have been proposed in the literature. 
But the process of combining two or more of these methods to form new ones have not been given serious 
research attention, therefore the motivation for this study. The aim of this study is to develop new structured 
methods of PI controller design for systems with time delay from existing methods. 

Combining existing methods of PI controller design to produce new methods will increase the 
number of methods in the PI controller design methods pool from which one can select. Also, depending on 
the system specification requirements, controller may be found in the new methods that will perform better 
than those in the existing methods. In view of these assumptions, the equations for the generation of stability 
region in the controller parameter space were derived. Using these equations and the critical frequency (we), 
the system stability boundaries were plotted. Based on the stability region, GA, WGC, and CCSR methods 
were used for designing PI controllers in the k,-k; convex stability region. The resulting controller gains 
from these methods were then combined in four different ways. Each of these ways stand for a new method 
and yielding four additional methods at the end. Making these methods depend solely on the results obtained 
from the existing methods. 

The existing and proposed methods were tested via simulation using unit step response. Considering 
the resulting system TDPMs value from the 3 examples used for demonstration, the proposed methods add to 
PI controller design methods pool. Also, some of the proposed methods outperformed some of the existing 
methods. 


2. METHOD 
2.1. Stability boundary for PI controlled system 

Figure 1 shows the block diagram of unity feedback control system used for the deriving of system 
transfer functions. The transfer functions for the plant, time delay and controller denoted as G,(s), Gq(s), 
and G(s), respectively are expressed in (1), (2) and (3). Based on these equations the system closed loop 
response (Y(s)) to the input (R(s)) can be determined. 


Figure 1. Block diagram of unit feedback time-delay control system 


Indonesian J Elec Eng & Comp Sci, Vol. 28, No. 3, December 2022: 1437-1450 


Indonesian J Elec Eng & Comp Sci ISSN: 2502-4752 O 1439 


_ Ns) 
Gp(s) = DO (1) 
Gals) =e (2) 


where the time delay in seconds is T, 


kps+ki 


G(s) = (3) 


sS 


where, kp and k; are the proportional and integral gains, respectively. 

Replacing “s” with “jw” in (1), and then to simplify the D-decomposition method, both the 
numerator N (s) and denominator D (s) of (1) were decomposed into their odd and even components as 
shown in (4). It should be noted that for compactness purpose (—w?) was removed from N,(—w?), 


N,(—@?), De(—w?) and Do (—w?) in (4). 


Ne+jwNo 
De+jwDo 


GpGw) = (4) 


The equations and the conditions required for determination of PI stabilization region in the kp-k; 
plane is presented as follows [26], [27]: 


For w = 0 

k; =0 (5) 

For w = 0 
_ (@*NoDot+NeDe) cos(wtT)+w(NoDe—NeDo) sin(wt) 

kp = —(Ne?+w2No*) (6) 
_ w*(NoDe—NeDo) cos(wt)—w(NeDe+w2NoDo) sin(wt) 

ki > —(Ne2+w2No*) (7) 

where w = [0, wc] 
The phase of the system of (1) when w = w, is 180°, this condition yielded (8). 
tan(wr) = 2NoPe“NePo) (8) 


NeDet+@2NoDo 


The value of w, is the smallest value of w at which the interception of the graphs of the left- and right-hand 
sides function of (8) against w occurred. The PI stability region in the kp-k; plane is the area bounded by the 
line produced by (5) for w = 0 and the stability locus generated by (6) and (7) for w = [0, we]. 


2.2. Determination of PI controller gains 

Three methods of PI controller gains determination which form the basis for the formulation of the 
new methods proposed in this study are considered in this section. These methods include GA, CCSR, and 
WGC which is presented in Subsections 2.2.1, 2.2.2 and 2.2.3 respectively. The objective of the GA is to 
optimize the PI controller gains (kp and k;). The objective of the CCSR is to determine the centroid of the 
stability region, while the objective of the WGC is to determine the weight of the geometrical center of the 
stability region. 


2.2.1. Determination of PI controller gains via GA 

The first step in the design of PI controller for a time-delay systems via GA is to formulate the 
optimisation constraints from the equations of the stability locus and the real roots line. Optimization of the 
PI controllers can be achieved by minimizing the objective function. In this study, integral of time multiplied 
by absolute error (ITAE) shown mathematically in (9) was used as objective function because it gives the 
best selectivity of the performance indices, it is a good criterion for PID controllers tunning [28] and when 
system parameters are subject to variations its minimum value is readily discernible [29]: 


ITAE = f,' tle(t)| dt (9) 


where e(t) is the error, t is the time in sec. and T, => T;,( T; is the system settling time in sec.). 
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Selection, crossover, and mutation which are the three reproduction operations used to form new 
generations in GA optimisation process are briefly explain as follows: 

a) Selection: Since the system stability region have been determined, the next step is to randomly select 
the initial population within the stability region. Evaluation of new generation is done by determining 
the objective function value for everyone in the population. The fitness function which was obtained by 
finding the reciprocal of the objective function, is normally used for ranking of the next generation in 
the population. Individuals are select for reproduction based on their fitness value. 

b) Crossover: This is the process by which part of the two parents in this case controller gains selected for 
the reproduction are combined to form new offspring to improve the next generations, 

c) Mutation: To guide against arrival at false optimum controller gains, mutation, which is the act of 
modifying individual gains in the generation, is applied. 

These three reproduction processes will continue until the adopted termination condition is satisfied. 

Detailed information on PID controller design for system with time delay using GA can be found in [21], [22]. 

Also, PID controller design for delay free systems is fully discussed in [30]. 


2.2.2. Determination of PI controller gains using CCSR method 

The coordinates of the corner and cusp points of the stability boundary locus were determined. 
Suppose the number of these corner and cusp points are n and m respectively, then their corresponding 
pj ki) and (kpg Kia) respectively, j=1,2,....n and 
q=1,2,....m.. The coordinates of the centroid of the convex stability region (kpcesr, Kiccsr) are the required 
controller gains in kp-k; plane, and they can be calculated by using (10) and (11) respectively [19]. 


coordinates in the k,,-k; plane can be represented by ( 


not m > 
L jar kp jt dias kpj 


kpccsr = (10) 


n+m 
1 ae 7 m > 
Lijat Kij+t}j=1 Kij 
n+m 


kiccsr = (11) 
2.2.3. Determination of PI controller gains using WGC method 

The generated convex stability region is made up of several points which can be defined by unique 
kp, k; coordinates. For m number of points on the stability locus, the k,, k; coordinates can be expressed as 
(Ky1, kis), (kp2, kiz), --- , (kpm, Kim). The stability boundary defined by k; = 0 is independent of œ and gives 
only real roots, therefore the coordinates of points on this line that corresponds to that on the locus are 
denoted by (Ky , 0), (Kp2, 0), ... , (Kpm, 0). Each of (Kp1, Kix), (Kp2, Kiz), --- , (Kpm> Kim) pair can be mapped 
to (kp1, 0), (kp2, 0), ... , (Kpm, 0), respectively. 

Note that the value of m is equal to the number of w steps (Nasp) and inversely proportional to step 
size (Aw). The expression relating 1,,5,, @, and Aw is shown in (12). Based on the stability boundary locus 
and real root boundaries mapping, the controller gains (Kpwgc and kiwgc) corresponding to WGC point of 
the system convex stability region can be determined using (13) and (14). This method is simple and can 
provide good transient performance for systems with time delay [23], [24]. Detailed discussion on WGC 
technique of PI controller design were reported in [23]-[25]. 


Nusp = 2 (12) 
dL. 

Kywec = op j= kpj (13) 
1 

kiwec = zp d= Kij (14) 


2.3. Proposed methods for PI controller design 

In the proposed technique, the controller gains were first obtained using GA, CCSR and WGC 
methods. Thereafter, the gains obtained from these three methods were now combined by finding the 
centroids of the lines joining any of the two gain locations, and the triangle whose vertices are the locations 
of the three gains within the convex stability region. Whichever the case either straight line or triangle, the 
resultant gains can be obtained by finding respective centroids in the convex stability region. 
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The four possible centroids in this case are: 
a) Centroid of the line joining the location of GA gains and CCSR gains, 
b) Centroid of the line joining the location of GA gains and WGC gains, 
c) Centroid of the line joining the location of CCSR gains and WGC gains and, 
d) Centroid of the triangle whose vertices are the location of GA gains, CCSR gains and WGC gains. 


2.4. Time domain-based performance analysis 

Performance evaluation of any designed control systems can easily be carried out via simulations. 
This is achieved in this study by plotting the step response of the controlled and uncontrolled systems on the 
same axes. Several system TDPMs have been defined for use in assessing system behaviours when the 
system test signal is a step function. The parameters used in this article include: i) rise time (T,,), 11) peak time 
(Tp), iii) percentage overshoot (%OS), iv) percentage undershoot (%US), v) settling time (T;) and vi) steady- 
state error (ess). 

Control system time domain response is made up of two states, namely the transient-state and 
steady-state responses. The transient-state response depends on system swiftness to the input signal and the 
degree of closeness of system response to the reference input. The former is characterized by T, and T, while 
the latter is dependent on T, and %OS. It should be noted that it is not possible to simultaneously satisfy same 
extreme condition for system swiftness to the input signal and closeness of system response to the reference 
input. Because these requirements are contradictory in nature, therefore compromise is normally reached to 
produce the best controller for a given system. The major parameter that is used for characterizing steady- 
state response is ess. 

In the actual sense, selection of controller for a given plant depends on the system performance 
requirements. If a system is required to reach the set point as fast as possible the controller gains should be 
selected such that the values of T, and T, are as lower as possible. But if tracking the reference input is what 
is important, then controller gains must be selected such that T, and %OS are as low as possible. It should 
also be noted that for good systems steady-state response, ess must be as low as possible. 


3. RESULTS AND DISCUSSIONS 

For comprehensive demonstration of the proposed methods, 3 different examples (first, second and 
third order systems) obtained from [23] were used. Since the steps involved in the PI controller design 
methods applied in this study are structured, detailed explanation on these steps were only present for 
Example 1. The detail is not present in Examples 2 and 3 to conserve energy and time. But detailed enough 
information that can be used for controller performance analysis and comparative analysis among the existing 
and proposed methds are presented. 


3.1. Example 1 

A first-order system with time delay was considered. The transfer functions of plant and delay are 
given in (15) and (16) respectively. Combining (15) and (16) with (3) will yield the system open-loop 
transfer function from which the closed-loop transfer function can be determined. Based on the closed-loop 
transfer function the system closed-loop characteristic equation can also be determined. 


Gp(s) = <5 (15) 
Ga(s) = e~ (16) 


From (2), (4), (15) and (16), t = Isec. Ne = 1, N, = 0, De = 1 and D, = 1. 

Substituting the values of tT. Ne, No, De, and D, into (8) and by letting the left- hand side and right- 
hand side functions be f,(w) and f (w) respectively, wç can be determined by plotting the two functions 
against w as shown in Figure 2. Note that w, is defined as the first point of intersection of the two plots, and 
from Figure 2, wç = 2.03 rad/sec. Using (5), (6), and (7), the system stability region of Figure 3 was obtained. 
It can be seen in Figure 3 that system convex stability region is bounded by the real roots line generated by 
(5) and the stability locus produced by (6) and (7). 

The determination of controller gains for Example 1 using GA was achieved by following the 
explanation of Subsection 2.2.1 The equation for the system stability locus is a polynomial of 6""-order and it 
was generated by curve fitting Figure 3, yielding (17) with R? = 0.9998: 


ki = pikp® + pak? + pakp* + paky® + pskp” + Pekp + p7 (17) 


New methods for proportional-integral controller design for time-delay system (Aye Taiwo Ajiboye) 


1442 O ISSN: 2502-4752 


where: p, = -0.030001, p, = 0.059705, p3p3 = 0.011301, p4p4 = -0.14149, psp5 = -0.26956, p,p6 = 0.93416, 
p7 = 1.1352. 

The constraints for GA were obtained from (5) and (17) while the objective function was based on 
the ITAE of (9). The GA was implemented, and the resulting PI controller, that is proportional and integral 
gains are kpga = 0.6649 and kjg4= 0.5519 respectively. As earlier explained in Subsection 2.2.2, in this 
example also, the controller gains at the centroid of system convex stability region were determined. From 
Figure 4, there are 2 corners and 1 cusp points on the stability boundary locus therefore, n=2 and m =1. 


2p 2r 
fl(w) 
1b 
1.57 
o fr(w) fl(w) 
=, N 
3 ~ ii ie 1 
z `~ adea 
S -1 F} Na ` z Stability region 
~ M 
2 > 0.5 F 
a 
2+ 7 
x203 / >< 
~ 
Y: -2.022 ES o 
~ 
-3f m- 
Tay < stablity locus 
eai — — ~ real root line 
-4 1 1 1 -0.5 1 i i 
o 0.5 1 1.5 2 2.5 3 3.5 4 S -0.5 0 0.5 1 15 2 25 
w (rad/sec.) kp 
Figure 2. Plots of fi(w) and f.(w) versus w for Figure 3. Plot of stability region for Example 1 
Example 1 
pae 
Stability region cusp points 
1.57 
1b 
z 
Convex stability region 
0.5 F 
Stability region corner points. 
0 
Stability locus 
— — — Real root line 
@ Corners and Cusp points 
05 1 1 1 1 A 1. i 
-1 -0.5 0 0.5 1 1.5 2 2.5 


Figure 4. Required points for determination of CCSR for Example 1 


The coordinate of the 2 corner points are (pis Ra) = (-1, 0) and (pas kis) = (2.263, 0) while that of 


the cusp point is (Can ka) = (1.1312, 1.717). Using (10) and (11) the coordinate of the centroid point is 
(kpccsr» Kiccsr) = (0.7981, 0.5723). WGC method of PI controller design is presented in Sub section 2.2.3. 
Considering the generated stability region of Figure 3 and the value of w,, Aw = 0.01 rad/s was selected. 
Using (12), Nysp= m=204. Substituting the value of m into (13) and (14) produces the controller gains 
Kpwec = 9.4442 and kiwgc = 0.4651 at the weighted geometrical centre point of the system stability region. 
The locations of (kpaa: Kiga ) = (0.6649, 0.5519), (kpccsr» Kiccsr) = (0.7981, 0.5723) and (kpwec, Kiwec) = 
(0.4442, 0.4651) are labelled A, B and C in the system stability region as shown in Figure 5, these locations 
form the vertices of triangle ABC. 
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Figure 5. Location of the GA, CCSR and WGC controller gains for Example 1 


Based on the coordinate of the vertices of the triangle ABC in Figure 5 and its side lines, the 


resultant gains at the four possible centroid points were determined by using the steps given in Subsection 
2.3. To find the centroid of a line, the number of corners will be 2 and the number of cusps will be O that 
means n =2 and m=0. Also, to find the centroid of triangle the number of corners will be 3 while the number 
of cusps will be 0 which implies that n =3 and m=0. 


a) 


b) 


c) 


d) 


The four possible centroid points for this example were determined using (10) and (11) as follows: 
Centroid of the line that joins the coordinate of GA gains (Kyga, Kica ) = (0.6649, 0.5519) and that of 
CCSR gains (kpcesr» Kiccsr) = (0.7981, 0.5723) that is line AB was determined to be (kpmo kimi) = 
(0.7315, 0.5621), 
Centroid of the line that joins the coordinate of CCSR gains (Kyccsr, Kiccsr) = (0.7981, 0.5723) and 
that of WGC gains (kKpwec, Kiwec) = (0.4442, 0.4651), line BC was determined to be (kpm2z Kim2) = 
(0.6211, 0.5187), 
Centroid of the line that joins the coordinate of GA gains (Ky¢a, Kiga ) = (0.6649, 0.5519) and that of 
WCC gains (Kpwec, Kiwec) = (0.4442, 0.4651), line CA was calculated and yielded (kpm3, kim3) = 
(0.5545, 0.5085) and 
Finally the centroid of the triangle with vertices defined by the coordinate of GA gains (Kyga4,Kiga ) = 
(0.6649, 0.5519), coordinate of CCSR gains (Kyccsr, Kiccsr) = (0.7981, 0.5723) and coordinate of WGC 
gains (Kpwec» Kiwec) = (0.4442, 0.4651) was also determined to be (kpm4 Kima) = (0.6357, 0.5298). 

The location of all the controller gains of the existing and proposed methods in the system convex 


stability region are shown in Figure 6. The performance of the proposed method was investigated via 
simulation. This was achieved by using unit step response of the system. The system was subjected to a unit 
step function after the incorporation of each of the 7 designed controllers (3 using existing methods and 4 
using the proposed methods). 
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point (kpWGC, kiWGC) “> 
point (kpCCSR, kiCCSR) 


e 
x 

a 

+ point (kpM1, kiM1) 7 
Pm point (kpM2, kiM2) 

Æ point (kpM3, kiM3) 
@ point (kpM4, kiM4) 


ki 
oO 
a 
N 
| 
| 
| 
= 
oO 
> 
fo) 


0.5 f 
0.48 F 
Cc 
x 
0.46 fi fi 1 fi f 
0.4 0.45 0.5 0.55 0.6 0.65 0.7 0.75 0.8 
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Figure 6. Location of all controller gains in the stability region for Example 1 
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3.2. Discussions of results of Example 1 

The determined controller gains using existing and proposed methods for example | are shown in 
Table 1. The corresponding unit step response for this system under the influence of the seven controller 
design methods and uncontrolled system (UCS) are shown in Figure 7. From Figure 7, the UCS response is 
far from the reference input, therefore the need for a controller. The required TDPMs for the analysis of the 
designed system were obtained from Figure 7. The parameters were used for producing the bar chart shown 
in Figure 8. 


Table 1. Existing and proposed methods controller gains for Example 1 
k,GA  k;GA k ,CCSR k;CCSR k, WGC k;WGC k M1 k;M1 k,M2 k;M2 k,M3 k;M3 k,M4 k;M4 
0.6649 0.5519 0.7981 0.5723 0.4442 0.4651 0.7315 0.5621 0.6211 0.5187 0.5545 0.5085 0.6357 0.5298 


tar 


a 
œ 


Response 
o 
[e>] 


0.4 


Time (seconds) 


Figure 7. Closed loop unit step response for Example 1 


From Figure 8, the minimum and maximum values of T, were obtained when CCSR and M3 were 
applied respectively. The minimum and maximum values of T, corresponds to CCSR and WGC respectively. 
Also, application of M2 and M1 yielded minimum and maximum values of T, respectively. But %OS was 
minimum when M2 was used and maximum for CCSR. ess = 0 for all methods of PI controller design used 
for this example. Based on these results, it can be inferred that CCSR gives the fastest response and M2 gives 
the best input tracking response. Since ess=o all the methods have excellent steady-state response. 


9 
OGA mwoec AccSR OM1 OM2 AM3 OM4 
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$5 A 
= A 
24 JA: 
a JAN 
E 4 
EAI 
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1 NA 
ar 
0 SEN 
Tr Tp Ts %0S ess 
TDPMs 


Figure 8. Bar chart for TDPMs and PI design methods for Example 1 


As can be clearly seen from Figure 8, compared to GA, M1 gives faster response while M2, M3, and 
M4 give better reference input tracking. Also, compared to WGC, M1 gives faster response and M2 give better 
reference input tracking. Lastly, compared to CCSR, methods M2, M3 and M4 produce better reference input 
tracking. 
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3.3. Example 2 
A second-order system with time delay. The plant and delay transfer functions are as expressed in 
(18) and (19) respectively: 


G(s) = Gas (18) 
Ga(s) = e* (19) 


where from (2) and (19) t = Isec. and from (4) and (18) N; = 1, N, = 0, De = 1 — w? and D, = 2. 

As done in Example 1, the location of all the controller gains in the stability region using the 
existing and proposed methods are shown in Figure 9. The performance of the proposed methods was 
simulated and investigated using system unit step response. The system was subjected to a unit step function 
after the incorporation of each of the 7 designed controllers (3 using existing method and 4 using the 
proposed method). 


ki 


— — —line AC 


05 0.55 06 065 07 075 08 085 09 095 1 
Figure 9. Location of all controller gains in the stability region for Example 2 


3.4. Discussions of results of Example 2 

The determined controller gains using the existing and proposed methods for Example 2 are shown 
in Table 2. In this example, the resulting unit step response for the seven controller scenarios understudied is 
shown in Figure 10. It can be seen from this figure that the response of UCS is far from matching the 
reference input which call for controller design. The time domain performance parameters associated with 
each of the controller design methods were obtained from Figure 10 and used for plotting the bar chart of 
Figure 11. 


Table 2. Existing and proposed methods controller gains for Example 2 
kpGA k;GA K CCSR Kk;CCSR k WGC k;WGC k, M1 k;M1 k, M2 k;M2 k, M3 kıM3 k,M4 k;M4 
0.6308 0.3427 0.9953 0.3803 0.5228 0.3087 0.8131 0.3615 0.7591 0.3445 0.5768 0.3257 0.7163 0.3439 


From Figure 11, the minimum and maximum values of T, correspond to CCSR and M1 
respectively. At the same time the minimum and maximum values of T, are associated with CCSR and WGC 
respectively. The minimum and maximum values of Ts were obtained using WGC and M4 while the 
minimum and maximum values of %0S were realized using WGC and CCSR respectively. The minimum 
ess was recorded when M1 was used while maximum value was obtained when either GA or M1 was used. 

The implication of these results is that the system of Example 2 will have the swiftest response if 
method CCSR is used but using method WGC will give a system with response that closely track the reference 
input. The system gives good steady-state response because highest percentage ess obtained was 1.28% which 
is less than the common standard of 2%. From Figure 11, it can infered that, compared to GA, methods M2, 
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M3, and M4 yield faster response, while M1 give better reference input tracking. More so, compared to WGC, 
methods M2, M3, and M4 gives fast response. Lastly, compared to CCSR, M1 gives better reference input 
tracking. 
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Figure 10. Closed loop unit step response for Example 2 
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Figure 11. Bar chart for TDPMs and PI design methods for Example 2 
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3.5. Example 3 
Design of PI controller for third-order system with time delay. The plant and delay transfer 
functions are shown in (20) and (21) respectively. 


1 
(s+1)(s+2)(s+3) 


G(s) = (20) 


Gals) = es (21) 


By comparing (2) and (21), tT = 0.5 sec. Also, comparing (4) and (20), Ne = 1, N, = 0, De = 6(1 — w°) and 
D, =11-0?. 

The locations of all the obtained controller gains in the stability region using the existing and 
proposed methods are shown in Figure 12. The performance of the proposed method was investigated via 
simulation using system unit step response. The system was subjected to a unit step function after the 
incorporation of each of the 7 designed controllers (the three existing and four proposed methods). 
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Figure 12. Location of all controller gains in the stability region for Example 3 


3.6. Discussions of results of Example 3 
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Table 3 shows the determined controller gains for the existing and proposed methods for Example 3. 
The resulting unit step response for Example 3 for the seven controller scenarios is shown in Figure 13. From 
this figure the response for the UCS is far from getting to the desired output, therefore the system needs a 
controller. The TDPMs needed for the analysis of this system were obtained from Figure 13. The parameters 


were used to plot the bar chart of Figure 14. 


Table 3. Existing and proposed methods controller gains for Example 3 


kpGA Kk;GA kpCCSR Kk;CCSR k WGC 


k,WGC 


kM1 k;M1 k, M2 k;M2 k M3 k;ıM3 


k,M4 


k;M4 


3.9360 2.6453 6.9330 3.1913 3.6823 


2.5755 


5.4345 2.9183 5.3076 2.8834 3.8091 2.6104 


4.8504 


2.8040 


From Figure 14 the minimum and maximum values of T, were obtained when CCSR and M1 was 
respectively used. Also, the minimum and maximum values of T, were associated with CCSR and WGC 
respectively. The minimum and maximum values T, were realized using WGC and CCSR, while the 
minimum and maximum values of %0S were produced using WGC and CCSR respectively. The value of 
ess was minimum when any out of GA, WGC and M1 was used but maximum when use is made of CCSR. 
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Figure 13. Unit step response for Example 3 
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Figure 14. Bar chart for TDPMs and PI design methods for Example 3 


The implication of these results is that the system of Example 3 will have the most swiftness 
response if CCSR method is used but using WGC method will result in a system with response most close to 
the reference input. Since the maximum percentage ess is 1.39% which is less than 2% which is the standard, 
the steady-state response achieved using any of the controller design methods is good. From Figure 14, when 
compared to CCSR, methods M2, M3, and M4 gives faster response, while M1 gives better reference input 
tracking. Furthermore, compared to WGC, methods M2, M3, and M4 give faster response. Finally, compared 
to CCSR, methods M1, M2, M3, and M4 gives good reference input tracking. 

In summary, any of the proposed methods (M1, M2, M3, and M4) can be selected, based on the 
system performance specifications. It should however be noted that it is not practicable to realize a system 
that will satisfy all the TDPMs requirements in real life due to their contradicting behavious. Therefore, 
tradeoffs are normally made once the system performance specification tolerance is known. 


4. CONCLUSIONS 

New PI controller design methods for systems with time delay were reported in this study. From the 
four structured PI controller design methods developed, the selection of any of the methods is informed by 
the system performance specifications. The system’s performance analysis was simplified based on its 
transient- and steady-state analysis. Satisfactory steady-state responses were obtained in each of the proposed 
controller design methods since the resulting percentage steady-state error is within the tolerable value (2%). 
The controlled systems were classified as swift and close reference input tracking according to the resulting 
system TDPMs. The performance of the controllers realised using few existing methods was tested against 
the proposed methods via simulation using closed loop unit step response. When compared to any one of the 
three existing methods, some of the proposed methods yield faster response speed at the expense of reference 
input tracking, while the reverse is the case for others. Any of the proposed methods (M1, M2, M3, and M4) 
is a suitable candidate depending on the required system performance specifications. One of the benefits of 
this study is the creation of pool for PI controller design methods from which selection can be made. 
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